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Abstract 

Fractals and multifractals and their associated scaling laws provide a quantification 
of the complexity of a variety of scale invariant complex systems. Here, we focus 
on lattice multifractals which exhibit complex exponents associated with observable 
log-periodicity. We perform detailed numerical analyses of lattice multifractals and 
explain the origin of three different scaling regions found in the moments. A novel 
numerical approach is proposed to extract the log- frequencies. In the non-lattice 
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sets of complex exponents of the lattice sequence converge to the set of complex 
exponents of the non-lattice one. An algorithm for the construction of the lattice 
sequence is proposed explicitly. 



* Corresponding author. Address: KPL F 38.2, Kreuzplatz 5, ETH Zurich, CH-8032 
Zurich, Switzerland. Phone: +41 44 632 89 17, Fax: +41 44 632 19 14. 

Email addresses: wxzhou@ecust.edu.cn (Wei-Xing Zhou), dsornette@ethz.ch 
(Didier Sornette). 

URL: http://www.er.ethz.ch/ (Didier Sornette). 
^ This work was partially supported by the NSFC (Grant 70501011) and the Fok 
Ying Tong Education Foundation (Grant 101086). 



Preprint submitted to Chaos, Solitons & Fractals 



2 February 2008 



1 Introduction 



The word fractal was coined by Mandelbrot to describe sets, which consist 
of parts similar to the whole and which can be described by a fractal dimen- 
sion [1], which is in most cases a fractional number. Meanwhile, multifractals 
were introduced to characterize the statistics of the strength of singularity of 
self-similar measures carried on general geometric supports [2|3|^f5ll6] . Frac- 
tals and multifractals have been extensively introduced in the analysis of a 
wide class of phenomena in physics, geophysics, chemical engineering, turbu- 
lence, growth phenomena, and so on. At the same time, new ingredients were 
added in the theory of multifractals. Based on the discovery of anomalies in 
DLA (Diffusion Limited Aggregation), the phenomenon of phase transition 
in the multifractal spectrum was found and studied [7f5|U] . This led to the 
mathematical framework of the so-called left-sided multifractals [TU]|11II12|[T3] . 
The vector-valued multifractals were introduced to formulate nonconservative 
measures such as the velocity field of turbulence [H]. Continuous multifractals 
p!5|16|17j were adopted to quantify the droplets breakup in turbulent jet flows 
|18] . which generalize the multipliers from discrete to continuous. 

Another important concept is the negative dimension or latent dimension 
p!9|20f21j . The physical implication of negative dimensions was discussed firstly 
in [22]. In practice, most of the physical processes are random, which leads to 
the phenomenon of sample-to-sample fluctuations of the multifractal function 
/(a) by an amount greater than what the error bars on any single sample 
would indicate. In general, negative dimensions arise in such random mul- 
tifractals [23|2l] . The randomness of multifractals is found at least in two 
cases where negative dimensions may emerge [23|2^ . First, one may obtain a 
multifractal constructed by a multiplicative cascade that is inherently proba- 
bilistic and the negative dimensions, if they exist, describe the rarely occurring 
events. Second, one may have to investigate the experiment from a probabilis- 
tic viewpoint. For instance, one-dimensional cuts of a deterministic measure 
carried out using a deterministic Sierpinski sponge will inevitably introduce 
randomness and may be regarded as random samples of a population. It is 
worth remarking that randomness in multinomial measures does not imply the 
existence of negative dimensions [23]. Actually, the multifractal slice theorem 
[26|27] presents a mathematical interpretation of negative dimensions and is 
the basis of experimental measurement of lower dimensional cuts when the 
measurement of the whole set is difficult to perform. The multifractal slice 
theorem makes rigorous some aspects of Mandelbrot's intuition that nega- 
tive dimensions may be explained geometrically by considering cuts of higher 
dimensional multifractals. 

In the context of fractals and multifractals, scale invariance is the most impor- 
tant concept of symmetry to characterize fractals or multifractal measures. In 
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its simplest form, scale invariance can be expressed by the following expression: 

0{r)=iO{Xr), (1) 

where O is the observable and A is the magnification factor or scaling ratio. 
In the common sense, A is a continuous parameter resulting in the so-called 
continuous scale invariance (CSI). In other words, CSI is associated with the 
fact that expression ([T]) holds for arbitrary values of A chosen in the set of 
positive real numbers. However, CSI is not always true when investigating 
fractals and multifractals. The general solution to Eq. (fl]) is not just a power 
law as CSI would imply but take the form 

0{r) = r-^^/j{r), (2) 

where D = ln£/lnA is the fractal dimension and ip{Xr) = ip{r) is a priori an 
arbitrary log-periodic function with period In A |28|29|3U|31II32|33] (the term 
"log-periodic" means that the function ip is periodic in the variable Inr with 
period In A). Actually, the recently proposed discrete scale invariance (DSI) 
and its associated complex dimensions are an elaboration of the concept of 
scale invariance in which that system is scale invariant only under powers of 
specific values of the magnification factor (see [51] for a review). DSI with the 
signature of log-periodic oscillations decorating power laws have been reported 
in several practical applications, such as in material rupture [35|36II37] , earth- 
quake precursors [38)139] . DLA [^0135] . turbulence [¥T|42|43] . and economics 
[44] . In addition, there are also reports in the theoretical aspects of fractals 
[28|29|3U|31|32|33] and multifractals [IS]. Notice the evolution of the concept 
of dimension from integer to fractional to negative and finally to complex 
numbers, each time with precise physical meaning. 

However, Eq. ([1]) is only a special (lattice) case of DSI (when A takes discrete 
values) [in]. In this work, we shall review the general framework of DSI in 
strictly self-similar fractals and multifractal measures in Appendix |X] with ex- 
tension to joint multifractal measures based on the self-similarity of moments. 
The structure of the complex exponents will also be discussed for both lattice 
and non-lattice cases. 

The rest of this paper is organized as follows. Extensive numerical simulations 
are carried out on lattice multifractal measures to investigate different scaling 
regimes stemming from finite size effects in Sec. |2l We propose a novel ap- 
proach for the extraction of log-periodicity in the lattice case in Sec. [31 Sec. jH 
concludes. The Appendix summarizes the general framework of discrete scale 
invariance in fractals and multifractals, which is used throughout the paper. 
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2 Numerical simulations on lattice multifractals 

In the absence of analytic results, this section presents computer simulations 
to explore the properties of lattice multifractal measures. The properties of 
the log-periodic function iplr) are investigated in details. This task remains 
difficult for lattice multifractals with more than three partitions in each gen- 
eration. With the increase of partition number n, the iteration number and 
hence the computational time needed in order to obtain detectable log-periodic 
structure with relatively low noise level increases exponentially. We thus have 
to confine our study to small n values. In Sec. 12. 1^ we study binomial measures 
with fixed amplification factors Ai = (-\/5 + l)/2 and A2 = (v^ + 3)/2 and 
varying measure multipliers mi = m and m2 = 1 — m to demonstrate three 
scaling regimes: the fractal regime, the crossover regime and the Euclidean 
regime, the latter resulting from finite size effects. Note that the geometric 
support is not fractal. Then, we study the properties of the Euclidean regime 
in Sec. 12.21 allowing us to gain insight into the structure of the noise stem- 
ming from the incomplete construction of the multifractal measure and of the 
partitioning. In Sec. 12. 3[ we study the scaling properties of the amplitude of 
log-periodic oscillations decorating the scaling of moments. 



2.1 Finite size effect 

We performed 13 iterations to construct a lattice binomial measure with mea- 
sure multipliers m = 0.2 and 1 — m = 0.8 (see [mfT3]l6] for the definition 
and construction of binomial and multinomial measures). Then the moments 
r(r, q) for different orders q were calculated based on the box-counting method 
with sizes r evenly distributed in logarithmic scale in [e~^, e~^] (see the Ap- 
pendix for definitions) . The log-log plot of F (r, q) against r for q = 1 recovers 
r(r, g = 1) = 1 (as it should from the normalization of the measure and the 
definition flA.9p of the moments) with minor fiuctuations of order 10~^^ stem- 
ming from the precision of the computation. The log- log plot of r(r, q) versus 
r for g = is a perfect straight line with slope —-Do = Recall that, if 
the geometric support is fractal, the slope should be — -D/. The log- log plot of 
r(r, q) against r for q G (0, 1) U (1, 00) looks like a strip with slope r(g) which 
is analogous to the situation presented in Fig. [TUJ We postpone the discussion 
of this case till Sec. O 

Fig. [T] shows the dependence of the moments F(r, q) with respect to box size r 
in a log-log plot for g < 0. There are three regimes clearly visible in the figure 
separated by two vertical dashed lines. Log-periodic oscillations decorate the 
scaling law in regime /. It is easy to calculate numerically the corresponding 
log-frequency /o ~ 1/lnA = 2.0781 by employing a spectral analysis on the 
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Fig. 1. The dependence of moments T(r,q) with respect to box size r in a log-log 
plot for q < 0. The multifractal is a lattice binomial measure with scale multipliers 
Ai = VA^= (V5 + 1)/2 and measure multipliers mi = 1 — 1712 = 0.2. We performed 
13 iterations to construct the multifractal. There are three regimes clearly visible 
separated by two vertical dashed lines, where / is the fractal regime, /// is the 
Euclidean regimes and II is the crossover regime. The left vertical line is located 



precisely at Tc 



,13 



'see text for details). 



residuals defined a.sip— (ip) or on the local derivatives defined in Eq. (1A.20P in 
the Appendix. Regime /// corresponds to perfect straight lines with different 
slopes. Regime // shows the crossover from large box sizes in regime / to small 
ones in regime III. We shall clarify in the sequel that / is the fractal regime, 
/// is the Euclidean regime and // is the crossover regime. We shall see that 
this phenomena is due to the finite size effect. 



It is well known that, although the mathematical model for the construction 
of a fractal or multifractal measure contains an infinite number of steps of 
the iterative construction process, the notion that they are (statistically) self- 
similar can only apply between certain cutoffs in practice [1]. In other words, 
what we are able to deal with in nature, laboratory and computer simulations 
are in fact pre-fractals and pre-multifractals. We nevertheless use the term 
"fractal" and "multifractal" without the prefix "pre" , following the custom in 
the literature. 



Let Tc denote the crossover scale (close to the inner cutoff) of a multifractal 
with measure /ij = /x(rj) distributed uniformly on a segment jFj of scale r^, 
where i = 1, ■ ■ ■ , n. It is clear that 

= (max{ri})\ (3) 



where b is the iteration number of the construction. This is the largest 
segment among the n'' segments with uniform measure. Choosing to count 
boxes with r < r^, we have approximately that r{F;r,q) ~ Y17=i^i^ij^yl)- 
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There are about Ui = Vi/r boxes covering segment J-'i, each box possessing the 
measure Hi/rii. Thus, we have r(jFj;r, g) = ni^Hi/riiy. It follows immediately 
that 

n 

r(F;r,g)^5:/i,r^V-\ (4) 

1=1 

where is independent of the box size r. We obtain that r(g) = q — 1 

from the definition flA.lSp and Dq = 1, corresponding to the Euclidian regime 
///. Concerning the general situation in a rf-dimensional space, we have r(g) = 
(g — l)d and Dq = dui regime ///. This result just states that a pre-fractal is 
Euclidean if we investigate it at a very high resolution. 

We have observed in Fig. [1] that Tc is precisely the right boundary of the 
Euclidean regime ///. The inner cutoff is obviously larger than r^, since a 
box with size greater than doesn't guarantee the absence of a possible 
destruction of the log-periodic structure close to r^. Indeed, such a box will 
inevitably split a segment into two parts which introduces errors. However, 
the inner cutoff can not be determined rigorously but only numerically. 

When changing the value of m, we can still observe four different kinds of 
moments r(r, g) associated with different g's: (1) g < 0, (2) g = 0, (3) g = 1 
and (4) g G (0,1) U (l,oo). It is apparent that the plots of the moments 
as a function of r for (2) and (3) are similar for all lattice and non-lattice 
multinomial measures, and thus trivial. For the present binomial measure, 
when m is large (close to 1), case (1) exhibits a strip while (2) shows three 
regimes. The results for small m on the three regimes also apply by symmetry 
m — i> 1 — m to large m — ^ 1. When m = (-\/5 — l)/2, the measure degenerates 
to a mono-fractal, where the measure is uniformly distributed over the support 
with density 1. 

2.2 Residuals of the moments in the Euclidean regime III 

Although regime /// is trivial in the sense that r(g) = (g — there are 
still fine structures stemming from those boxes covering two segments and the 
rightmost box as well. Let us define residuals in regime /// as 

0(r, g) = In r(r, g) — (g — 1) Inr. (5) 

For different g < 0, the residuals 0(r, g) have similar shapes. Figure [2] plots 
the residual for g = —20 with respect to r. We performed 6 = 5 iterations to 
construct the binomial measure. We see accelerating oscillations with decaying 
amplitudes towards r = 0. The maximal box size is = 0.0902. We labelled 
this point as 0. The consecutive points of local minima and local maxima are 
labelled as 1, 2, ■ ■ ■ in turn as shown in Fig. [2l The coordinates of these points 
are denoted accordingly (r^, ^(ri, g)) with i = 0,1,2,---. The n-th oscillation 
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Fig. 2. The dependence of the residual 0(r, q = —20) defined by expression ([5]) with 
respect to r. The even numbers are used to label the points of local maxima, while 
the odd numbers denote points of local minima. The point is precisely at r = rc- 
It is visible that 4>{r2n,q) is constant for a given q. 

is from point 2n — 2 through point 2n — 1 to point 2n. One can see that the 
oscillations are similar to each other. We find that 

0(r2n,g) ~ 0(rc,g). (6) 

Note that Tq = Vc- The function (f){rc, q) is found to be linear with respect to 
q with negative slope. 

Figure [3] draws the box sizes r2n-i and r2„ versus n + 1 in a log-log plot. We 
find the following power laws: 

r2n-i = 2ri(n + l)"i (7) 

for n G and 

r2n = rc{n + iy^ (8) 

for n G Z~. It is important to note that r2„_i and r2„ are independent of q for 
given binomial measures. For different iterations b, and ri are not constant 
anymore. We can determine r2„ analytically by combining Eqs. ([3]) and ([8]), 
while Ti should be calculated numerically for different b. Let (5r)„ = r2„ — 
r2„_2, then (5r)„ = l/n(n + l). When n — > oo, it follows that (5r)„+i/(5r)„ 
1, which means that the oscillations are approximately periodic locally for 
large n. This is not surprised since (5r)„, decays more and more slowly when 
r ^ following a law. 

Define the amplitude of the n-th oscillation as the distance of point (r2n-i, 4>{'f2n 
to the line 0(r, g) = (0(r2„,g))„: 

hn{q) = 0(»^2n-l, q) - {(P{r2n, q))n, (9) 
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Fig. 3. The power law relationship of r2n (with n G Z~) and ?'2n-i (with n G Z'^) 
with respect to n + 1. The scaling exponents are both —1. 
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Fig. 4. The power law dependence of the amplitude hn{q) with q = —20 versus r2n- 
The scaling exponent is a = 1.0569. 

where ()„ is an average over different oscillations n. Figure H] shows the de- 
pendence of the amplitude hn{q) for q = —20 with respect to r2n- Again, we 
observe a remarkable power law relationship: 

hn{q) = ciq) ■ r^^, « = 1.0569, (10) 

where 

c{q) = h{q){rj2)-^ (11) 

is a function of q. 

We also investigated the dependence of hi{q) with respect to q, as shown in 
Fig. O It seems that hi{q) approaches a constant when q tends to — oo. The 
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Fig. 5. The dependence of hi{q) with respect to q. There is a limit for hi{q) when q 
tends to — oo. The inset shows the power law dependence of /ii(— oo) — hi{q) versus 
—q. The scaling exponent is /3 = 1.0013. 

points (o) in Fig. [5] can be represented by the following law: 

h{q) = h{-oo)-C/\qf, (12) 

where the values oo) = 0.6915, j3 = 1.0013, and C = 0.6293 are ob- 
tained from a nonlinear regression procedure. The inset verifies the power 
law dependence of oo) — hi{q) as a function of —q. We also find that 
oo) = 0.6915 is universal for different b and m, while /3 and C vary for 
different b and m. In addition, the amplitude hi{q) increases with 



2.3 Residuals of the moments in the fractal regime I 



We now study the log-periodic function ip{r, q) and the corresponding residuals 
ln'?/'(r, q) = lnF(r, q) — T{q) Inr of moments in regime / defined in Eq. f lA.17p 
of the Appendix. Denoting Ai = A, we find A2 = A^. From (lA.34p . we have 
7729_\'r(9) + (1 — niY X'^'^'^''^ = 1 which is a solvable univariate second-order equa- 
tion with respect to A"^^"^-*. The explicit expression of the exponent r(g) follows: 



In 2 - In 



m"? + x/m^? + 4(1 - m) 

tx 



where A = (v^ + l)/2 is the preferred scaling ratio. The behavior of ip is 
simpler, since it is strictly periodic in Inr. Thus, given m and q, the log- 
periodic oscillations have constant amplitude as a function of r. Figure [6] shows 
the dependence of ip{r,q) for q = —20 with respect to Inr/ In A for different 
measure multipliers m. 
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Fig. 6. The dependence of 'ijj{r, q) for q = —20 with respect to Inr/ln A for different 
measure multiphers m. The periodicity is clearly visible with period In A in the Inr 
variable. The amplitude of the oscillations decreases with increasing m for fixed q. 
The local minima are located precisely at r = A~* with i £ Z~ . 

The three lines in Fig. [6] are periodic in Inr / In A with period 1. The amplitudes 
of the oscillations decrease with increasing m for fixed q. The local minima are 
located precisely at r = A~* with i G Z~ . The abscissa of the local maxima are 
different for different m, while those for same m but different q overlap. For 
large r, the noise increases with m. For small r, the oscillations are gradually 
spoiled to the right of the crossover regime //, which were not shown in Fig. [61 
Hence, it is better to investigate the oscillations in the middle of the fractal 
regime / to obtain a quantitative description of the log-periodic oscillations 
of ip{r, q). 

As shown in Fig. El the local minima minr'?/^(r, q) are close to 1 at r = A~* with 
i G . Actually, in the trivial case r = 1, the moments r(r, g) are always 
equal to 1 for all kinds of conservative multifractal measures and all orders 
q. We calculated the values of min^ '?/'(r, g) for different q and m. The results 
are shown in Fig. [71 We find that min^ '4){r, q) approaches closer and closer to 
1 with the increase of |g| or the decrease of m. The distance of the minima 
miii.^ 'ipir^q) to 1 is exponential as a function of q for fixed m. The slopes of 
the straight lines fitted to the data for each m in Fig. [71 are approximately 
evenly spaced for small m. This verifies the analytic results that = 1. 

Figure [HI shows the local maxima ma.Xj.ip{r,q) versus q for different m. Note 
that moments of larger m contain a large noise amplitude making the numer- 
ical determination of local maxima rather difficult. We find that maxr^(r, q) 
increases exponentially with This result still holds for other lattice multi- 
fractal measures. The slopes of the straight lines fitted to the data for each m 
in the log-linear plot of Fig. [HI are proportional to Inm. 
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Fig. 7. The exponential approach of min^ •ip{r, q) to 1 when \q\ increases for different 
m is shown. The slopes of the straight hnes are approximately evenly spaced for 
small m. 
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Fig. 8. The exponential dependence of max^ ^l^{r, q) as a function of q for different 
m. The slopes of the straight lines in this log-linear plot are proportional to In m. 

We define the amplitude of log-periodic oscillations as 



h{q) = maxip^r, q) — mmip{r, q). (14) 

It is easy to see that h{q) has a behavior similar to that of ma:Xril^{r,q). It 
follows the important result that the amplitude of the log-periodic oscillations 
of lattice multinomial measures increases with A similar result is known 
for (lattice) fractals [30] • This may be related in a qualitative way to the ob- 
servation that large log-periodic oscillations are found in DLA and in rupture 
processes compared to very small log-periodic amplitudes in spin models de- 
fined on hierarchical geometries. Indeed, fractal growth and rupture processes 
are controlled by local field / stress concentration, i.e., by the high-order mo- 
ments of the field involved in the construction of the patterns. 
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3 A novel approach for extracting log-frequency 

A natural question arises in practice concerning the determination of the 
log-frequency and thus the preferred scaling ratio in the detection of log- 
periodicity in fractals or multifractals. Up to now, five different methods have 
attempted to improve the detection of log-periodicity [13]. The common part 
of four non-parametric methods among the five methods is the Lomb peri- 
odogram analysis. The Lomb analysis of ^/'(r, g) in Fig. [H] gives a very signif- 
icant peak at log-frequency /o = 1/lnA. Here, we study more complicated 
measures for which the log-periodic oscillations are not visible to the bare 
eyes but can still be detected by spectral analysis. We develop in this section 
a novel method for the detection of log-frequencies, which is expected to be 
powerful also more generally in the search of log-periodicity in real systems. 

We construct a triadic multifractal measure which has three partitions in each 
generation with scaling ratios Ai = 2 and A2 = A3 = 4 and measure multipliers 
nil = ^3 = 0.5 and m2 = 0. The measure is supported by the generalized 
Cantor set. Thus the preferred scaling ratio is A = 2 and the fundamental log- 
periodic frequency is /o = l/ln2. The analytic value of r(g) is also given by 
(1131) with m = 0.5 and A = 2. We use n = 11 iterations for the construction of 
the measure. Again, the counting box size is evenly sampled in the logarithm 
of the scale. The moments r(r, q) of the measure were thus calculated for 
different orders q based on the box-counting method. 

We present in Fig. [9] a typical plot of lnr(r, q) versus Inr for g < 0. Although 
the plot is very noisy, a linear least square fitting gives a slope close to the 
analytic exponent T{q). No unambiguous simple (log-)periodic oscillations are 
visible. Since Inrc = — nlnAi = —7.62, we restricted the interval of study to 
—6.17 < Inr, where the lower bound is indicated by the left downward arrow 
in Fig. O Indeed, a behavior different from large r can be observed for small 
r. We also introduced an outer cutoff in the large r range as indicated by the 
right upward arrow in Fig. [9l Our analysis for the detection of log-periodicity 
was thus carried out in the interval confined within the two arrows, containing 
a total of A = 131 data points. Note that a different but reasonable choice of 
this interval doesn't impact on the results. 

Fig. [To] shows a typical plot of lnr(r, g) versus Inr for g > 1. A linear least 
square fit of the points recovers a numerical estimate of r(g) which is close to 
the theoretical value. The plot also exhibits many parallel patterns (rectangle- 
like) which are arranged periodically in Inr, as indicated by the arrows. We 
find that the period is approximately 2 In A, which is interpreted as the fre- 
quency of a subharmonic to the fundamental log-frequency. A similar sub- 
harmonic log-periodic frequency was also reportedly detected in the growing 
diffusion-limited aggregation clusters [10]. However, the Lomb periodogram of 
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Inr 

Fig. 9. The dependence of lnr(r, q) with respect to Inr for q = —20. Although the 
plot is very noisy, a linear least square fitting gives a slope close to the analytic 
exponent T{q). No unambiguous (log-)periodic oscillations are visible. 

the detrended lnr(r, g) in Fig. [10] doesn't provide convincing signal of such 
/o/2 subharmonic (not shown). We nevertheless observed in the Lomb peri- 
odogram two outstanding peaks at log-frequencies / = 11.54 and / = 15.87, 
although this observation does not ensure the statistical significance of log- 
periodicity. Lower order moments provide more significant Lomb periodogram 
when q > 0, which means that the noise was amplified rapidly for large g's. 
On the other hand, the impact of q on the noise is less important when g < 0. 
This noise amplification effect of larger q suggests to investigate lower-order 
moments to detect log-periodicity. 
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Fig. 10. The dependence of lnT{r,q) with respect to Inr for q = 20. The arrows 
indicate a periodic pattern in Inr. The period is 2 In A which is interpreted as a 
subharmonic of the fundamental log-frequency. 

We detrended the moment lnr(r, g = —3) to remove the power-law scaling. 
Here the scaling exponent is estimated directly from the data points and the 
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Fig. 11. The Lomb periodogram of the residuals of the q = —3 order moment. 
The left inset shows the local details around / = 9.5/ In A, which is the symmetric 
axis of the periodogram. The right inset shows the local details around / = 19/ In A, 
beyond which the Lomb periodogram duplicates itself periodically. The Lomb peaks 
indicated by open circles are approximately evenly spaced. 

resultant residuals have zero mean. This method was used extensively in prac- 
tice (this method is different from that used in Sec. 12.21 and 12. 3p . Analyzing 
these residuals by the Lomb method gives a Lomb periodogram shown in 
Fig. [TTJ It is interesting to point out that the Lomb periodogram is periodic 
of period 19/ In A = 27.4112. The local detail around / = 27.4112 is shown in 
the right inset of Fig. [TTl We see a clear burst exactly at / = 27.4112 causing 
the discontinuity at this point. It is necessary to stress that this burst itself is 
not a Lomb peak. It is thus sufficient to analyze only the / < 27.4112 part of 
the periodogram as shown in Fig. [TTl In addition, the Lomb periodogram in 
Fig. [TTjis apparently symmetric with respect to / = 9.5/lnA = 13.7056. The 
left inset shows the local detail around / = 13.7056, which is a discontinuous 
burst as well. This is the reason why the right inset is symmetric with respect 
to / = 27.4112. 

In Fig. [TT|, the two highest peaks are located at / = 11.54 and / = 15.87. This 
is what we can see in the Lomb periodogram for g > 1 as we alluded above. The 
claim that the fundamental logarithm frequency is /o = 11.54 or /o = 15.87, as 
done in many other practical systems, has a false-alarm probability of 0.006% 
assuming that the noise is independent and Gaussian The existence 

of two peaks of the same height and the symmetric structure of the Lomb 
periodogram reinforce further the statistical significance of the existence of a 
genuine log-periodicity. 

An additional important evidence is offered by expanding ^/^(r, q) in the Fourier 
space, which identifies a discrete set of log-frequencies fn with amplitudes 
An. Here, we find that the amplitudes An of different frequencies are not 
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Fig. 12. The linear dependence of log-frequencies of the harmonics of the funda- 
mental frequency /o as a function of their order n. The slope of the fitted straight 
line is 1.4396 it 0.0005 providing an estimate of /o in excellent agreement with the 
theoretic value of /o = 1/ In 2. The preferred scaling ratio is thus A = 2.0030±0.0005 
having only a 0.3% deviation from the true value of A = 2. 

restricted to a presumed inequality that Ax > An with n > 1 (we refer to 
[50] for a detailed discussion). In the case where several An are comparable, as 
shown in Fig.lTT], we have to study the harmonics of the fundamental frequency. 
There are indeed many slim peaks in Fig. [H] that are approximately evenly 
spaced. We indicated these 18 peaks with open circles, whose corresponding 
log-frequencies are denoted as /„. This recognition process is objective, since 
all except the 6th and 13th peaks are higher than the remaining peaks. It 
is quite reasonable to interpret /„ as the n-th harmonic of the fundamental 
log-frequency. A convenient estimate of /o is the mean of the differences of 
consecutive frequencies. This gives /o = (/« — fn-i) = 1.4457 ib 0.0671. In 
addition, we can expect that /„ ~ nfo- Thus, we can take the slope of the 
line fitted to points (n, /„) as another estimate of /q. Figure [12] plots /„ as a 
function of n, which exhibits an almost perfect linear relationship. The slope 
of the fitted line is /o = 1.4396 ± 0.0005. The estimates of /o from the two 
approaches are comparable to each other and both in excellent agreement 
with the theoretic value of /o = 1.4427. However, the linear regression method 
of Figure [12] is obviously much more stable. The preferred scaling ratio is 
thus calculated to be A = e^^^ = 2.0030 with a standard deviation ax = 
aje^'f IP = 0.0005. 

In real systems, it is possible that the Lomb periodogram is much more noisy 
than found in Fig. [TT] In that case, we may only determine a few but not all 
harmonics, so that the successively detected harmonics are not consecutive. 
In this case, taking the differences between successively detected harmonics 
will not give a constant value equal to the fundamental log-frequency. This 
certainly increases the difficulty of the detection of log-periodicity and low- 
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ers the significance of tlie extracted log-frequency. One should thus estimate 
subjectively the possible order of the detected harmonics and reconstruct the 
best match according to a hnear fit such as in Figure [T^l but with holes for the 
missing harmonics. Notwithstanding its difficulties, this method was success- 
fully applied to the search of log-periodic in the energy dissipation of three 
dimensional fully developed turbulence [13]. 



4 Concluding remarks 

We have reviewed and generalized in this paper the basic discrete scale in- 
variance equations concerning the moments r(r, q) of fractals, multifractal 
measures and joint multifractal measures based on the presence of (statisti- 
cal) self-similarity in these objects. There exist a set T of complex dimensions 
for fractals and complex exponents for multifractals, corresponding to the set 
of complex solutions of the corresponding DSl equation (DSl: discrete scale 
invariance). The set T lies in a strip [r'^(g), ro(g)] x TZ in the complex plane 
and is symmetric with respect to the real axis, where T^{q) G T is the only 
exponent on the real axis. 

In the lattice case, the complex exponents are located evenly with spacing u; 
(where u is the fundamental angular log-frequency) on finitely many vertical 
lines in the complex plane. Corresponding to the complex exponents on the 
vertical line r(g) = r^(g), there is a log-periodic expression of r(r, g) with a 
fundamental log-frequency / = 1/lnA, where A is the preferred scaling ratio. 
We have verified the log-periodicity in multifractals numerically and proposed 
a novel precise approach for the extraction of /. Our simulations also show that 
the signal-to-noise ratio of the moments increases with |g| when g < or g > 0. 
Since larger signal-to-noise ratio corresponds to more significant Lomb peak, 
we should investigate lower-order moments to extract log-frequencies, which 
is independent of the fact that higher-order moments have larger amplitudes 
of oscillations. 

In the non-lattice case, the complex exponents are irregularly located in the 
strip and T^{q) is the unique complex exponent in the vertical line r(g) = 
T^{q). There are thus no preferred scaling ratio and no exact log-periodicity. 
However, we can construct a sequence of lattice multifractals whose preferred 
scaling ratios converge to 1 from above to approximate the non-lattice mul- 
tifractal so that the set of these lattice multifractals converge to the set of 
complex exponents of the non-lattice multifractal. In the Appendix, we pro- 
pose an explicit algorithm for the construction of the lattice sequence. 

Therefore, we can recognize lattice fractals and/or multifractals in real systems 
from the evidence of the presence of log-periodicity. We can obtain a preferred 
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scaling ratio A based on spectral analysis. This search not only sheds new light 
towards a better understanding of the system under study but also offers a 
constraint on the modelling of the system. Of course, A does not contain a 
detailed information on all the scaling ratios Aj. It should be noted that the 
search for log-periodicity has been carried out already in several real systems, 
such as in the growing diffusion- limited aggregation (DLA) clusters [ID], in 
two dimensional free decaying turbulence [12] and in three dimensional fully 
developed turbulence [5T|l3] . It is interesting that the preferred scaling ratio 
is A ~ 2 for the cases of DLA and 3D turbulence and it is possible that the 
log-frequency /o = 4 ~ 5 in the two dimensional free decaying turbulence case 
[12] turns out to be a harmonic of /o = l/ln2. However, we have to stress 
that the existence of a preferred scaling ratio A does not imply that there is 
only one scaling ratio equal to the measured A. For instance, A = 2 found in 
3D turbulence does not mean that the p-model [S2] is the correct and only one 
model. We shall further clarify this point by taking the growing DLA clusters 
as an example as follows. 

There are many models proposed to explain the fractal growth of DLA clus- 
ters, such as the hierarchical ghost model [53|5l] and the Fibonacci model 
[5^56|57j . In the hierarchical ghost model, pairs of single particles collide and 
aggregate to form dimers. Then pairs of dimers collide and aggregate to form 
four-particle clusters. This hierarchical growth process goes on and on to form 
eventually a large DLA cluster. Defining b{N) as the average number of bonds 
along the connecting pathway between two particles on an iV-particle cluster, 
we obtain the DSI equation [53] 

6(2iV) = ^b{N). (15) 



Since the average radius i? of a A^-particle cluster obeys R ~ yb{N) and 
~ R^, we have D = 21n2/ln(3/2) and the set of complex exponents with 
respect to (ITS]) is 

T = {2/D + 2nm/ ln2, ne Z}. (16) 

It is clear that A = 2 and the set of complex exponents is on a vertical line 
with abscissa of 2/D in the complex plane. On the other hand, by applying the 
results of the wavelet transform modulus maxima (WTMM) representation to 
the DLA azimuthal Cantor set, it follows that the number N{a) of maxima of 
the WTMM skeleton at scale a satisfies the Fibonacci rule 



N{a) = N{Xa) + N{X'^a). (17) 
It is easy to show that the set of complex dimensions is [^0|l6] 



„ r Inch 2?™^ r Ind) (2r2 + l)7r^ 
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where n E Z and (j) = (-\/5 + l)/2 is the golden mean. We see that there are 
two vertical lines in the complex plane in this model. Numerical simulations 
[inj lead to the identification of two log- frequencies /o = 1/lnA and fo/2 in 
excellent agreement with the prediction of the Fibonacci model. The preferred 
scaling ratio was estimated numerically to be A = 2.2. It shows that the 
Fibonacci model with two scaling ratios A and A^ provides a better description 
than the hierarchical ghost model in the sense of discrete scale invariance. 
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A General framework of discrete scale invariance in fractals and 
multifractals 

In this Appendix, we revisit the equations governing the concepts of discrete 
scale invariance (DSI) and the associated complex exponents in fractals and 
multifractal measures, and generalize them to joint multifractal measures. We 
shall show that the set of complex exponents lie symmetrically with respect 
to the real axis in a vertical strip in the complex plane. Furthermore, in the 
lattice case, the set of complex exponents lie on finitely many vertical lines 
periodically and the DSI leads to the log-periodic corrections to scaling im- 
plying a preferred scaling ratio that is proved to be independent of the order 
of the moments. In the non-lattice case, there are no visible log-periodicity, 
i.e., no preferred scaling ratio since the set of complex exponents are spread 
irregularly within the complex plane. A non-lattice multifractal can be ap- 
proximated by a sequence of lattice multifractals so that the sets of complex 
exponents of the lattice sequence converge to the set of complex exponents of 
the non-lattice one. An algorithm for the construction of the lattice sequence 
is proposed explicitly. 

A . 1 Revisiting DSI in fractals 

Smith et al [30] studied three Cantor sets with the initiators being 101 (the 
classic Cantor set), 101001001, and 101010001. They followed the Grassberger- 
Procaccia (G-P) approach [58] and investigated the periodic oscillations in 
the log- log plot of the correlation integral C{€) with respect to scale In 
this section, we revisited DSI in the classic Cantor set using an alternative 
approach [28] . 

In order to generalize to other self-similar fractal sets, consider a fractal set JF, 
which is the union of n disjoint parts Ti , namely T = VS^^iTi with Ti^Tj = $ 
for i 7^ j. Let N{T\r) be the number of boxes covering T with scale r. We 
assume that the magnification of Ti by a factor Aj results in T. This leads to 



Note that the reciprocal of the amplification factor, 1/Aj, is called the scale 
multiplier. Since 





n 





i=l 



we obtain 



n 




(A.3) 



i=l 
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Equation (1A.3P has a solution if and only if there exist integers ki and A 
satisfying 

\ = A^% (A.4) 

which is the generalization of the Fibonacci rule of the number of maxima of 
the WTMM skeleton at a given scale in diffusion-limited aggregation |3D] and 
corresponds to the so-called lattice case [16] . The solution can be expressed as 



N{T;r) = r-'^ip{r). (A.5) 

where '0(r) = ?/'(Ar) and D is the fractal dimension determined by the real 
solution of equation 

n 

EA-^ = 1. (A.6) 

i=l 

We see that N{T] r) exhibits log-periodic oscillations with respect to r. Since 
N{J^, 1) = 1, it follows that ?/'(A~^) = 1, where £ G Z~. For = 1, we recover 
Eq. O. 



A. 2 DSI m multinomial measures 



A. 2.1 General equation of DSI in multinomial measure 

We construct a self-similar multinomial measure ^ supported by a set JF by 
a repeated n-based multiplicative cascade, where can be an Euclidean set 
or a fractal set of dimension dimjF. In the first step, the support T carrying 
some measure is partitioned into n small pieces jFj with i = 1, ■ ■ • , n based 
on some given set of scale amplification factors 

A.. =11^11/11-^. II, (A.7) 

where the measure is redistributed on J^i according to a set of given multipliers: 

mi = n{J^i;r/\i)/n{J^;r). (A. 8) 

In the next step, each piece jFj is further partitioned into n smaller pieces 
according to the same amplification factors Aj and its carrying measure is 
redistributed again with the same multipliers mj. This procedure continues ad 
infinity. So far, we have constructed a deterministic self-similar multinomial 
measure. Define a g-order moment function 

r(^;r,g)=^/i'^(^;r), (A.9) 

r 

where /U.(JF; r) is the measure in an arbitrary box of scale r and J2r represents 
the sum over all boxes which result from partitioning with scale r. When 
we use boxes of size r to cover J^, the minimum number of disjoint boxes is 
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[II II /''"]• O'^^ finds easily that there might exist a box covering a part of T 
whose size is less than r. In this case, we can arrange the [|| || /r] boxes in a 
row (like a string) and cover T in infinitely different ways by translating the 
boxes string along JF, thus introducing a random "phase." This phase is one 
source of noise in the box-counting method. When r is small enough, the noise 
is reduced greatly so that the moments corresponding to different phases are 
approximately identical. In other words, the moment function is not affected 
by the partition procedure and we have 

n 

r(^;r,g)=^r(^,;r,g) . (A.IO) 



Substitution of Eq. (lAlSD in Eq. ([Qj) leads to 



r(^,;r,g)=mfr(^;V,g), (A.ll) 

which shows the self- similarity of the moments. Combination of Eqs. f lA.lOp 
and (lA.lip gives 

n 

Y{T;r,q)=Y.Y[T;\r,q)ml (A.12) 

This is the basic DSI equation of the deterministic self-similar multinomial 
measures. For statistically self-similar measures [59] , Eq. (IA.12P can be revised 

as 

fe n 

V[T\ r, g) = X! H r(JF; A,jr, q)mljPj, (A. 13) 

j=i i=i 

where pj is the probability of choosing the j-th rule out of the k rules and Aj j 
and accordingly the scaling ratios and measure multipliers [25] . 



A. 2. 2 Solution of the DSI equation for lattice multifractals 

In the case of g = 0, the moment function equals to the number of partitioned 
boxes r(jF; r, 0) = N{J^] r), which recovers Eq. flA.Sp . 

If Ai = • ■ ■ = A„ = A, it follows from Eq. f[02|) that 

n 

r(^;r,g)=r(^;Ar,g)^m?. (A.14) 

i=l 

In the formalism of multifractals, one expects that 

r(J^;r,g) oc r"(«), (A.15) 
where the scaling exponent of the moment is in the form 

n 

r(g) = -InJ^m^/lnA. (A.16) 
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In general, the solution to Eq. (1A.14P is 

r(J^^;r,g) = r'^'VW, (A.17) 

where ^l:{r) = ip{Xr) and r(g) is defined in Eq. (IA.16p . Since T{T] 1, g) = 1, 
we have = 1, where £ G Z~. 

Now, consider the case of lattice multifractals with the scale multipliers Aj 
satisfying condition flA.4p . It follows from flA.12p that 

n 

r(^;r,g)=^r(^;A'='r,g)m?. (A.18) 
1=1 

It is easy to verify that Eq. (]A.17p is also the general solution to Eq. (lA.lSp . 
where r(g) is the unique real solution of 

n 

Y^X'^'-^i^m'! = 1, (A. 19) 
1=1 

We emphasize that ip is related to q but independent of A. Also, the gen- 
eral solution of Eq. (lA.lSP has the same form of (1A.17P with different "T(g)- 
generation function" in the lattice case. 



A. 2. 3 Oscillations of other characteristic functions of multifractals 

It is natural to investigate the presence of oscillations for other functions 
characterizing multifractals, for instance, the generalized dimensions Dq{r), 
the strengths of singularity a{r,q), and the multifractal function f{r,a) or 
/(r, q) [3]|^ll5]|6] . when one investigates scaling laws in some experimental data 
as well as simulated data. Our analysis is performed in terms of the variable r 
since all these functions vary with r. As usual, we define the exponent function 
of the moments as 

T{r,q) = d\nT{r,q)/d\nr. (A.20) 
In equation (lA.20p . we remove the symbol of the support and use rir.q) 



instead of r(g) as is used in [31H] . Simple algebraic calculations on equation 
( lATH) leads to 

r(r,g) = r(g) + ^(r,g), (A.21) 

where 



Similarly, we have 



D,M^lim4^ = fl,+ limHM), ,A.23) 

qi^q q' — I q'^q q' — I 
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l) = — ^ = + 

and 

/(r, g) = qa{r, q) - T{r, q)f{q) + q — ^(r, q). (A.25) 

Hence, the four characteristic functions of multifractals, Dq{r, q), T{r, q), a(r, q), 
and /(r, q), are all log-periodic in r with the same period In A, which is inde- 
pendent of the order q. This result was verified in an experimental situation 
by extracting the log-periodicity of the moments of the energy dissipation in 
three-dimensional fully developed turbulence [13] . We shall return to this issue 
in Sec. IA.4I 



A. 3 DSI in joint multifractal measures 



In complex systems, there may be more than one self-similar measure embed- 
ded in the same sample space or phase space. For instance, in heated turbulent 
flow, both the energy dissipation rate e and scalars, say concentration C or 
temperature T, are known to be multifractal. Thus, it is interesting to con- 
sider the joint multifractal measures. Suppose that /ij with j = 1, ■ ■ ■ , Z are I 
self-similar measures embedded in T. The joint moments of order q = I]j=i Qj 
of the / measures are given by 



J{J^;r,q) = Y.]lf^fi^;r). (A.26) 

r j=l 



The measure multipliers of /ij are defined by 

rriij = Hj{J^;r/\i)/nj{J^;r). (A.27) 
Combining flA.26P and flA.27p . we have 

J(^,; r, q) = J(^; V, q) [] m^^j- (A.28) 

i=i 

It follows from the additive rule of the joint moment 

n 

J{T;r,q)=Y.J{T,;r,q) (A.29) 

i=l 

that the DSI equation is 

n I 

J{T; r, g) = ^ J{T; A^r, q) J] m:i!j, (A.30) 
i=i j=i 
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for deterministic measures where Xij and rriij are the i-th multiphers of fij or 

n' n I 

J{T; r,q) = J2Y^ .J{T; A,,,.r, q) [] mg.fcPfc, (A.31) 

k=l 1=1 j=l 

for statistically self-similar measures where Xt^k and Triij^k are i-th multipliers 
of fij for A;-th rule with probability p^. 

Similarly, if fij satisfy the lattice condition so that Xi^k = A'^'-'', the general 
solution of Eq. (lA.3ip is 

J{J^;r,q) = r'^^^'- ''?')^/;(r), (A.32) 

where ip{r) = ip{Xr) and r(gi, ■ ■ ■ , g^) is the unique real solution of the follow- 
ing equation: 

n' n I 
k=l i=l j=l 

Again, since J(JF; l,g) = 1, we have ip{X~^) = 1, where i G Z~. The deter- 
ministic case Eq. ( ]A.30I) can be viewed as a special case of Eq. (]A.3ip with 
n' = 1. 

It is clear that the DSI equations of multifractal measures in Sec. IA.2I are 
special cases of those of the joint multifractal measures with / = 1. Moreover, 
their DSI equations are the same in essence despite the possible difference in 
the coefficients. We expect to confirm the relevance of these formulas in real 
systems with hierarchical structures, especially in turbulence. 



A. 4 Complex exponents 



The previous sections emphasized the log-periodic structure of lattice fractals 
and multifractals. However, there is no log-periodicity for the non-lattice case. 
Instead, we can adopt the more general concept of complex exponents as 
was done in the former literature pn]|34f46] . Indeed, log-periodicity is only a 
special case of complex exponents in the lattice case, while complex exponents 
exist for both lattice and non-lattice fractals and multifractals. The complex 
dimensions of fractal strings were extensively studied from a mathematical 
view point [60)1^ . We note that a fractal string is not identical to the original 
fractal but to its complement. Hence, we propose a different approach, even 
for the fractals discussed in the previous sections. In this section, we study 
the structure of complex exponents in the complex plane and then provide an 
explicit algorithm for the construction of a sequence of lattice multifractals 
approximating a non-lattice multifractal whose set of complex exponents is 
the limit set of the sets of complex exponents of the lattice sequence. 
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A. 4-1 The structure of complex exponents 



In the previous sections, we assume that r(g) is real, thus the prefactor of the 
scahng relation r(r, q) oc r'^^'^^ is a function of r in the lattice expressed 
in Eq. flA.171) . Expanding ip in the Fourier space yields complex exponents |34j . 
Therefore, we can assume that r(g) is complex so that the prefactor of the 
scaling relation is independent of r. In this sense, we have 

n 

Em?A[(^) = l. (A.34) 

i=l 



A more general formulae is Eq. (]A.33p . The solutions of Eq. (1A.34P in the 
complex plane are denoted as 

r{q) =rR{q)+iTj{q), (A.35) 



where i = ^/—l is the imaginary unit. We shall study this Eq. (]A.34p in the 
sequel. The results are the same for the stochastic case. 

For convenience, let /i(r) = Y17=i Since oo) = oo and h{oo) = —1, 

and because /i(r) is continuous, real solutions to (]A.34p exist. Moreover, /i(r) 
is differentiable and we have dh{T)/dT < 0. Thus Eq. (lA.34p has a unique real 
solution, denoted as T^{q) hereafter. The geometric properties of r^(g) are 
well-known [6]. As for the complex solutions, it is easy to show from the real 
value of /i(r) that, if r(g) is a solution to flA.34p . then r(g) = TR^q) — «r/(g) 
is also a solution. Therefore the set of complex exponents is symmetric with 
respect to the real axis (the complex exponents come in pairs of complex 
conjugates). 

Let r(g) be a solution of Eq. (]A.34p . Then r^i^q) > T^{q). Otherwise, by 
assuming that TR{q) < T^{q), we have | ELi ^[^"^1 < Y.7=i X^''^'''^ < 
Xf^'^ = 1. This means that the complex exponents lie on the right 
of the vertical line r^^q) = T^{q). In the lattice case, r(g) = T^{q) + inu 
also belong to the set of complex exponents, where n & Z and uj = 27r/lnA 
is the fundamental angular log-frequency of lattice fractals and multifractals. 
In the non-lattice case, r(g) = r^(g) is the only complex exponent satisfying 
'''nil) = '^^{q)- Furthermore, there is a right boundary ro(g) of the set of 
complex exponents, that is, TR^q) < ro(g). (We refer the interested reader to 
Sec. 2.5 of |l6] for a rigorous proof for fractal strings, which can be easily 
generalized to the present case.) Hence the set of complex exponents lie in a 
strip [T^{q),To{q)] x 7^. 

In the case of lattice multifractal, there exists a preferred scaling ratio A so 
that equation (]A.34p can be expressed as a polynomial equation (]A.19P with 
respect to the unknown A^'-^-' of degree maxj{kj}. The set of complex exponents 
is obtained easily from the set of complex solutions of (1A.19I) . Hence there are 
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finite complex numbers Tj{q) such that the set of complex exponents is given 

by 

r = {Tj{q)+inuj, neZ, J = l,2,...,i}, (A.36) 

where £ < maxj{kj}. Thus the complex exponents lie on finitely many vertical 
lines in the complex plane. The points on each vertical line corresponding to 
a given Tj{q) are distributed evenly with spacing u. 

Assume that (tj)/? 7^ {Tk)R for j k and Ti(g) = T^{q). Then we have 

{t,)j = uj/2, (A.37) 

where j > 1. This result derives because the points on each line are evenly 
located with spacing u and symmetric with respect to the real axis. The log- 
frequencies of the moments are thus given by 

;^^„ = Mi_±i!^i=( = \ (A,38) 

[{l + n)u/2n, if J > 1 

where n & Z — {0} when j = 1 and n ^ Z when j > 1, f\^\ is the fundamental 
log-frequency and /j>i,o are subharmonics of f\^\ [ID]. Certainly, it is also 
possible that there are at least two complex exponents with indices j and k 
such that {jj)R = {Tk)R. In this case, the points on the line TR^q) = {rkjR are 
no longer evenly spaced with period uj. This possibility will be clearer in the 
following when we approximate a non-lattice multifractal with a sequence of 
lattice multifractals. It is clear that /j>i,o might change for different orders q. 
Nevertheless, /i^i = uj/2tt holds for any q. In other words, the fundamental 
log- frequency of a lattice multifractal is independent of the order q [13] . 

In the case of non-lattice multifractals, the set of complex exponents is not 
periodic. However, the necessary and sufficient condition for the existence 
of log-periodicity is the existence of periodicity in the plane of the complex 
exponents associated with a preferred scaling ratio. In the present case, there 
is no such preferred scaling ratio and hence no log-periodicity. The complex 
exponents of a non-lattice multifractal can be approximated by the complex 
exponents of a sequence of lattice multifractals with larger and larger uj. The 
approximation sequence of lattice multifractals can be constructed explicitly, 
as in the case of fractal strings [16]. This construction process is described in 
the following Sec. IA.4.21 



A. 4-2 Approximating non-lattice multifractal by a sequence of lattice multi- 
fractals 

Let us denote the scaling ratios A of the non-lattice multifractal as where 
j = I,-- - ,n. Given any real number Aq > 1, we have a set of numbers 
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6 = {6j = InAj/lnAo : j = 1, ■ ■ ■ ,n}, defining the coordinates of a point in 
tlie n dimensional space 7?.". Since tlie multifractal is non-lattice, there is at 
least one j such that 6j is an irrational number for any possible choice of Aq. It 
is known that there exists a sequence of points e^'^) = {ef^ : J = 1, . . . , G 
7^" which converge to 6, where Oj''^ are rational numbers. Obviously, the 
sequence of rational numbers 6'^-^^ converge to Oj for every j = 1, - ■ ■ ,n. If 
some 6j is rational, we stipulate trivially that = 6j for all k. Posing 



Aj'^^ = Aq^ , we obtain a sequence of lattice multifractals A'^'^^ with scaling 



g(fe) 



ratios {Aj'^^ : j = 1, ■ ■ ■ ,n} converging to the non-lattice multifractal. The set 
of complex exponents T^''^ of A^^^^ converges to the set of complex exponents 
T of A. In practice, we can take a lattice multifractal A^'^^ with k large enough 
as an approximation of the non-lattice multifractal A. 

Let us write = /t^p where and t^p are relatively prime. Let t^^^ 

be the minimal common divisor of {t^!^^ : j = 1, ■ ■ ■ Then A^^-* = Aq''*' ' is 
the preferred scaling ratio of A'^'^^ We thus recover the polynomial equation 
dAUD: 

1 '7 

= 1, (A.39) 



Ao 



where sfk^'^^tf^ are integers. We can choose a sequence ensuring that 6) 
is monotonous with respect to k. It follows that t*^'^-* increases with k. There- 
fore, the fundamental log-frequencies /'•'^'^ = l/lnA'-'^^ of A^^^ increase when 
converging to A. In other words, there is no finite limit of Z*^^-* when k tends 
to oo and the fundamental log-frequency of a non-lattice multifractal is thus 
not defined. Practically, when one investigates the Lomb periodogram of the 
residuals of a g-order moment of a non-lattice multifractal, the log-frequency 
corresponding to the highest Lomb peak is indefinite and increases when the 
sampling of data points of the moment becomes denser (so as to investigate 
more scales r) [13]. On the other hand, the distance of separation between 
the complex exponents on each vertical line becomes larger and larger. This 
is not surprising since the degree maxjjSj'^^t'^'^V^i'^^} '^^ polynomial equa- 
tion (1A.39P increases with k implying more vertical lines, some of which may 
overlap. 
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